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E. B. Frost, Director of the Yerkes Observatory, Williams Bay, 
Wisconsin, U.S.A.; 

The Rev. J. G. Hagen, S.J., Director of the Vatican Observa¬ 
tory, Rome; and 

Johannes Franz Hartmann, Astrophysical Observatory,-Potsdam, 
Germany. 


Seventy-one presents were announced as having been received 
since the last meeting, including, amongst others:— 

J. G. Bohn, die Kunst-Uhren auf der K.-K. Sternwarte zu Prag, 
presented by the Prague Observatory; The History of the 
Geological Society of London, by H. B. Woodward, presented by 
the Society. 

4 Sixteen charts of the Astrographic Chart of the heavens, 
presented by the Royal Observatory, Greenwich. 


On the Lunar Inequalities due to the Motion of the Ecliptic and the 
^ Figure of the Earth. By Ernest W. Brown, Sc.D., F.R.S. 

1 1 . The general disturbing function for the motion of the 
ecliptic. —Let O l9 0 2 , 0 3 be the angular velocities of a set of 
moving rectangular axes about themselves; x,y, z the co-ordinates; 
u , w the velocities of a particle with respect to these axes; and 
let F be the force function divided by the mass of the particle. 
Then the equations of motion are given by 

du n , /) 0F 

dv , j /j 9F ' 

Tt - w e l+ ue^--, 


where 


Put 


where 


dw /j , /, 0F 
~dt - u6 2 + v6 1 = — , 


v =^- s e 1+x e s , 

w= %~ x0 ‘ i+y6 ^ ’ 


H = \{u 2 + v 2 + w 2 ) - F - R , 

R = tx9 z - wx0 2 + wy6 1 - uy0 3 + uz0 2 - vzO l ; 
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then, if we assume that 9^ 0 2 , 0 3 are independent of x , y, z , u, v , 
w, we may write the equations of motion in the canonical form, 


(ho 

- 0H 

dx_d H 

dt 

dx 9 

dt du 5 

dv ___ 

0H 

^/ = 0H 

dt 

dy’ 

dt dv ’ 

dw __ 

0 H 

dz_d H 

dt 

02 ’ 

dt dw 


If E be neglected, the equations become of the same form as 
for fixed axes; hence R is the disturbing function for the motions 
of the axes. These motions are so small that we can neglect 
squares and higher products of the terms in R, and therefore can 
treat R as a disturbing function of the same nature as that used 
for ordinary planetary perturbations. 

In order to obtain 9 V 0 2 , 0 3 , we put 1 for the inclination of the 
moving ecliptic (xy plane) to the fixed ecliptic (that of 1850 * 0 ), 
r for the longitude of the node on the fixed ecliptic, L for the 
angle which the zx plane makes with the plane containing the 
poles of the fixed and moving ecliptics. Then, by Euler’s 
equations, 

/1 d% • f • ./ t dT 

0 = — sin L-sin^ cos L—-, 

1 dt dt 

/I dt r .; • -r C?T 

— — cos L + sin ^ sm LP- , 

2 dt dt . 


^3 


dr , 

= cos l + 
dt 


<dL 

dt 


As L is at our disposal, we so take it that the distances of the 
origins of reckoning ori the fixed and moving ecliptics from their 
-common node are the same, that is, so that L = 90 ° - t. The 
moving axis of x therefore passes through a “ departure point.” 

When we are given the values of i, r in terms of the time, we 
have the material necessary for the solution of the problem. 


2 . As a matter of fact, ^-are so small that their squares 

dt dt 

and products may be neglected. We can therefore put 


sin i’ *= 1 , cos / =i ; t = const., 


and obtain 


/) dt 

e ^di C0ST ’ 


/) dt' . 

0 2 = -smT, 


0 3 = o. 
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Whence 


E = { (wy - vz) cos r + (uz - wx) sin r |, 

dx • 

or, since u differs from by a quantity of the same order as i , etc., 

at 

tj di' {( dz dy\ * ( dz dx\ . ) 

^ = At\\ y dt- Z I) C 0 ST -\ X dt- Z dt) SmT ) 

di' 

= suppose. 
at 


The method of the variation of arbitrary constants permits us 
to substitute in Q the values of x , y f z , and of their derivatives in 
terms of the arbitrary constants and the time as found in the 
solution of the main problem of the lunar theory. These values 
can be obtained directly from the results which I have given in. my 
papers in the Memoirs of the Society on this subject. 

3. The value of ( is of the form 

i' = pt + P, 

where p is a constant and P a sum of periodic terms whose 
arguments depend on those of the motions of the Earth and 
planets. The latter terms have generally been neglected, but I 
shall show that they give rise to a few terms which are not 
insensible. At this stage the question arises as to what moving 
ecliptic we shall refer the motion of the Moon. We may refer it to 
the actual ecliptic, in which case I have found by calculation that 
P gives rise to a number of primary terms of short period and a 
few of long period. Or we may refer it to the mean ecliptic, i.e. 
put P = o in the above formulae; but if we do so, it will bo 
necessary to introduce the latitude of the Sun above this mean 
ecliptic into F. With this latter method, I have found that all 
the primary short-period terms practically disappear and the long- 
period terms have rather smaller coefficients than with the former 
plane of reference. 

I therefore adopt the mean ecliptic together with the one or two- 
very minute terms which are of long period relative to that of the 
Moon’s node, and which do not then give rise to any terms in the 
Moon’s co-ordinates. Hence, for the computation, we may put 

E 1= i5Q. 

4. The principal part of the force function for the Moon’s, 
motion under the influence of .the Earth and Sun is 

2+z2) 1 

r r s [ 2 r 2 2 ) 

where x\ y\ z are the co-ordinates of the Sun. 
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The main problem is solved with z = o. Hence, the additional 
portion due to z is, if we neglect squares of z\ 

-p __m' 3 (xx' + yy')zz' 

-^2-' 

With the notation and limitations of § 2, 3 we have 


z = i'(y f cos r — x sin r). 

d 2 x 0 F 

Also, with this value for z and by means of the equations —= = — 

dt 2 dx j 

etc., it is easy to show that 


Hence, since t is a constant, we obtain 


E»2 — 


.,dQ 

dt 


i 


where for 1 we substitute its periodic portions, omitting those 
mentioned at the end of § 3. 

5. It is of some interest to compare the two disturbing 
functions R 1? R 2 . In general we have 


Thus the effect of the rotation of the axes alone, apart from the 
motion of the Sun out of the fixed plane, is to introduce a term 




From these expressions we may also prove the state¬ 


ments made in § 3 with reference to the presence of long- and 
short-period terms, remembering the two forms which have been 

given for ^ 5 .. In fact, the terms of long period in Q are relatively 


much larger than those in ^ owing to the presence of small 

lit { 

divisors ; for a term of very long period in R 2 , the portion arising 

from ~(i' Q) is small compared with that arising from R,. 
dt 


6. The computation from the two functions R p R 2 , as defined 
in §§ 3, 4, give the results which follow. The values of Leverrier * 
have been used for i\ t, and the method is the same as that which 
I have used for planetary inequalities in general. 

I put w v io 2 , w s for the mean longitudes of the Moon, of its 
perigee and of its node; Z, D, F, n , e, y are the quantities a£ 
defined by Delaunay; T, Y, J are the mean longitudes of the 


* Ann. Obs. Paris (Afdm.), vol. iv. pp. 13-21, 
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Earth, Venus, and Jupiter, respectively; cf) — w 8 + 96°'2; ft — w x 
+ 96°*2. 

Then the variations of the elements for the principal inequality, 
argument <£, arising from R 1} are 

8w 1 = - o"*289 sin <^>, +0"*001471 cos <£, 

8w 2 —4- o"*840 sin <£, 8e = *ooo cos <£, 

8 w b = - 1 s '’59 sin 8 y= + o"*698 cos <£ . 

These agree very closely with the values found by Professor 
Newcomb.* They differ a little more from those found by 
Dr. Hill; t he, however, used a method with a literal develop¬ 
ment, and the small differences are probably due to slow con¬ 
vergence. 

• The principal term in latitude, argument <£', has a coefficient of 

+ i"*382. 

7. The new terms arising from R 2 are, if we omit terms with 
coefficients less than o"*oio, 

In Longitude. In Latitude. 

+ o"*oi9 sin (<£ + 5T -3V+119 0 ) -fo"*o77 sin (<£' + 5T- 3V + n9°‘4) 

+ o"*oo3 sin (<£ + 2V) 2 - 2J 4 * 90°), + o"*o30 sin - 5T + 3V - H9°‘4) 

+ 0#, ‘° 3 S sin (<t>' + 2J + 72 0 ) 

+ o"*oi8 sill (<£' - 2J - 72 0 ). 

The two terms in longitude are “ primary ” terms of long 
period. The second one is of interest, .since its period is that of 
the empirical term which appears to represent best the long-period 
difference between theory and observation, namely, about 280 
years. Its coefficient is, however, too small to account for this 
difference. J There is a term with the same argument due to 
the direct effect of the planets, but its coefficient is also insensible. 
The terms in latitude are “ secondary,” and have periods approxi¬ 
mating to a month. 

8. Inequalities arising from the figure of the Earth .—These 
have been computed very fully by Dr. Hill,§ who used the method 
of Delaunay, and there would be no reason for the mention of 
them here if it were not that slight doubts have arisen concerning 
their degrees of accuracy on account of slow convergence. I have 
used the equations of variations and the results of my lunar theory 
which allow the coefficients to be obtained easily within o"*oo3. 
This method is brief, and the computations occupied about eight 

* Carnegie Institute, publ. 72, p. 132. The difference in the value of tin 
is due to a difference in the definition of this variation. 

t Annals of Math., vol. i. p. 57. Coll. Works, vol. ii. p. 77. 

t In a note on p. 170 of the present volume, and in this paper as read 
before the Society, I gave o"*2i as the value of this coefficient. A recalcula¬ 
tion with the first form of R2 revealed an error in the former computation. 

§ Amer. Eph. Papers, vol. iii. pp. 201-344. Coll. Works, vol. ii. pp. 
181-320. 
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days. I put i j/ = w 8 -b the precession, and obtain for tbe variations 
of the elements 

Sw 1= + 7"*317 sini/r, Sn = - o"*oo9 cos \[r , 

8 w 2 = — 2 ,f 'og2 sin ifr , 8 e — + o"‘oo2 cos ij /, 

8 w 3 — + <) 6 "' 6 () sin if/ , 8 y = - 4"*351 cos l/s . 

The shorhperiod terms add - o"*oi7 sin (\f/ + E) to the latitude. 

The principal term in longitude is and the principal elliptic 
terms are + o' A, 519 sin (t // + l) + o"*515 sin (i/r - T). The principal 
term in latitude is - 8"*355 sin + E), and that with argument 
ijf- F is + o"* 338 sin fj/ - F). When Hill’s results are reduced 
to the value of the ellipticity used above (1/296*3),* my coefficient 
of the principal term in latitude agrees with his, and my coefficient 
of the principal term in longitude is o' 7, 030 less. The coefficients 
for the term with argument ij/ - F also agree; this fact furnishes a 
useful test, since this coefficient is the difference of two numbers 
each nearly thirteen times .as large as the coefficient. The two 
principal elliptic terms are o"*02i, o^oiy, respectively, greater 
than those of Hill. An examination of Hill’s literal developments 
for the coefficients shows that the small differences can all be 
explained by slow convergence. 

The complete results for the classes of terms considered here 
will be given in the fifth (concluding) part of my Theory of the 
Motion of the Moon , with the terms due to the action of the 
planets and to perturbations not considered in the previous parts. 

New Haven t Conn.: 

1908 February I. 


Postscript (1908 April 4).—The terms of the second order have 
been examined and produce nothing.sensible. However, in conse¬ 
quence of my attention being called by Dr. Hill to a doubt as to 
whether there was a portion of the secular acceleration due to the 
figure of the Earth terms being referred to a moving ecliptic, I made 
an actual computation of this second order perturbation and found 
that the greatest effect could be exhibited in' the form of a term of 
period about 15,000 years and coefficient o"‘i5. This, equivalent 
to a secular acceleration of o"*oooi within historic times, is entirely 
insensible. The computation will be given in chapter xiv/of my 
memoir just mentioned. 

* Corresponding to the result marked (0) in M.N ., vol. lxiv. p. 531. 
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